Flat-band ferromagnetism proposed for an organic polymer crystal 
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Motivated from the flat-band ferromagnetism conceived theoretically for a single chain of five 
membered rings (polyaminotriazole) by Arita et a/., [Phys. Rev. Lett. 88, 127202 (2002)], we 
have studied whether the magnetism can indeed occur as a bulk, i.e., in a three-dimensional crystal 
of the polymer, by means of the spin density functional calculation. We find that the intra- chain 
ferromagnetism is robust against crystallization as far as the flat band is made half-filled. We have 
further investigated the actual crystal doped with various compounds, where HF2 is shown to put 
the system close to the bulk ferromagnetism while stronger anions such as BF4 or PF6 should be 
promising. 
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PACS numbers: 75.10.Lp, 71.20.Rv, 71.10Fd 

Materials design, especially that from the standpoint 
of many-body effect, is one of the most challenging av- 
enues in condensed matter physics. In particular, de- 
signing magnets in materials consisting entirely of non- 
magnetic elements is an issue of great interest. In fact, a 
variety of ideas for exotic ferromagnetism have been pro- 
posed in this decade, where first-principles calculations 
are combined with electron-correlation studies for models 
exemplified by the Hubbard model [1-5]. In this context, 
organic ferromagnets are of special interest. While or- 
ganic ferromagnetism has been realized in localized-spin 
systems such as TDAE-C60 [6], whether we can have a 
band (i.e., itinerant) ferromagnetism in or ganics remains 
a challenging question. 

Recently, we have proposed a new possibility of itin- 
erant ferromagnetism in an organic polymer of five- 
membered rings, polyaminotriazole (PAT), and have ex- 
plained the origin of the magnetism in terms of Mielke 
and Tasaki's flat-band ferromagnetism [7] by coupling a 
first-principles band calculation with a study for the Hub- 
bard model [4]. Concisely, the flat- band ferromagnetism 
occurs when there is a dispersionless band which satisfies 
a special condition (called connectivity condition) in the 
one-electron band structure. Note that the magnetism 
is distinct from the usual narrow-band limit in that the 
magnetism occurs as a quantum interference effect in the 
repulsively interacting system for a finite transfer energy. 

While the monomer (aminotriazole) is commercially 
available, experimental attempts at its polymerization 
have just begun [8]. Experimentally, we have to aim at 
a bulk magnetization, so theoretically it is imperative to 
study whether we can have a ferromagnetism in three- 
dimensional crystals, which is of fundamental as well as 
practical interests. 

The key issues are (i) the robustness of the intra- 
chain ferromagnetism to start with (i.e., whether the 
inter-chain interaction in the crystal disturbs the intra- 



chain ferromagnetism), (ii) whether the inter-chain in- 
teraction can be ferromagnetic (i.e., the competition be- 
tween the ferromagnetic (direct-exchange) inter-chain in- 
teraction and the inter-chain energy gained by the an- 
tiferromagnetic order), and their energy scale, and (iii) 
which chemial dopants should be opted to bring the flat- 
band close to half-filling, a necessary condition for the 
fiat-band ferromagnetism. The present paper addresses 
these problems. 

Problem (i) is non-trivial, since the presence of adja- 
cent chains can affect the magnetism within the chain [9] . 
Conversely, if the inter-chain interaction does not destroy 
the intra-chain ferromagnetism, this is already interest- 
ing for experimental opportunities, since if we can for 
example epitaxially grow the organic crystal on a metal- 
lic substrate with a large work function such as platinum, 
the interface state of the PAT chain can possibly have a 
net polarization. In fact, we shall show that the presence 
of neighboring chains in the crystal does not destroy the 
intra-chain ferromagnetism for the realistic chain-chain 
distance. 

On the other hand, problem (ii) on the inter-chain 
magnetic coupling, we shall show that ferromagnetic and 
antiferromagnetic interactions are both small in magni- 
tude and their competition is subtle. 

As for (iii), we we have studied various anions as 
dopants in PAT crystal to examine the magnetic prop- 
erty. We shall conclude that, while the dopant studied 
here, CIO4, F, and HF2, bring the system close to, but 
not in, the ferromagnetic phase, the chemical trend indi- 
cates that stronger anions having higher electron affinity 
(BF4 and PFe) are expected to realize the flat-band fer- 
romagnetism. 

We have adopted a first-principles band calculation 
within the framework of the generalized gradient ap- 
proximation based on the spin-density functional the- 
ory (GGA-SDFT) to compare the total energy of vari- 
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ous magnetic states. In the GGA-SDFT calculations the 
exchange-correlation functional introduced by Perdew, 
Burke, and Wang [10] is adopted with the ultra-soft pseu- 
dopotential [11,12] in a separable form. The wave func- 
tions are expanded by plane waves up to a cut-off energy 
of 20.25 Ry. 

To clarify the physics the doping has been done in 
two ways. We first look at the band- filling dependence 
by reducing the number of electrons, where a uniform 
negative background charge is introduced to make the 
system charge- neutral as in Ref. [4]. The latter part of 
the paper examines actual anions as dopants, where the 
charge transfer across the anion and PAT is calculated 
self-consistently. 

We start with the atomic configuration of the 
(undoped) PAT crystal obtained by the structure- 
optimization, shown in Fig.l, where the total energy is 
minimized with the conjugate gradient scheme [13]. We 
have assumed that the unit-cell is orthorhombic and con- 
tains two chains of PAT which form a herringbone struc- 
ture. We have studied two types of the crystal struc- 
ture, A(non-staggered) and B (staggered). Similar crys- 
tal structures have been considered in recent calculations 
for polythiophene [14] and poly-phenylene-vinylene [9], 
where PT2 and PT1 in Ref. [14] correspond to our A 
and B, respectively. 




FIG. 1. The optimized atomic configuration of crystallized 
polyaminotriazole (before doping) for type A and type B. The 
green, white, blue balls represent C, H and N, respectively, 
and the square represents the unit cell. 



number of electrons. We focus on the case of the half- 
filled flat band, for which the flat-band ferromagnetism is 
originally conceived. In the present crystal, top four va- 
lence bands are nearly flat to be precise (corresponding to 
the presence of four five-membered rings in a unit cell), 
and the half-filling refers to the case of one carrier per 
ring. We have obtained four types of solution for both A 
and B crystal structures in the doped case: an intra-chain 
ferromagnetic (F) and inter-chain antiferromagnetic(AF) 
state (which we call inter AF hereafter), an intra-chain 
AF and inter-chain F state (intraAF), a state antiferro- 
magnetic both for intra-chain and inter-chain (AF 2 ), and 
finally the ferromagnetic state. 

The band structure for the doped system is shown in 
Fig. 2 for the ferromagnetic solution. For both A and B 
structures, we can see that the dispersion along the chain 
is quite similar to that of an isolated PAT chain (shown 
in the inset), from which we can see that the effect of 
the inter-chain interaction is small for the ferromagnetic 
state. 

hole doping (structure A) 
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FIG. 2. The band structure for the ferromagnetic solution 
for the hole-doped crystallized polyaminotriazole in type A (a) 
and type B(b) structure. The black (green) lines represent the 
bands for the majority (minority) spin. Top-right inset shows 
the band structure for an isolated chain of doped polyamino- 
triazole. 



We have determined the size of the unit-cell as fol- 
lows. We first determined the dimension of the unit-cell 
along the chain to be 7.17 A, which minimizes the ground 
state energy of an isolated chain [4] . For the unit-cell size 
perpendicular to the chain, we have calculated the total 
energies by changing the linear dimension with the in- 
crement of 1.0 au (i.e., 0.529 A) to search for the energy 
minimum. The resulting unit-cell size is 8.46Ax 6.35A. 

We have then doped the system, first by reducing the 



The total energies of the three AF states for the struc- 
ture A(B), as measured from those of the F state in re- 
spective structures, are, respectively : £?(interAF) = —40 
(0) meV, ^(intraAF) = 130 (100) meV and E(AF 2 ) = 40 
(100) meV. If we estimate the "intra-chain magnetic cou- 
pling" as Jintra = (intraAF) — i£(F), and "inter-chain 
magnetic coupling" as Jinter = i^(interAF) — E(F), Jintra 
is ~ —100 meV (< 0, i.e., ferromagnetic) for both A and 
B structures. If we compare this value with that for a sin- 



gle PAT chain (estimated to be ~ -50 x 2 = -100 meV 
[4], where the factor 2 is required for comparison with 
the present value for two chains per unit cell), the intra- 
chain ferromagnetism is seen to be quite robust even in 
the crystal. 

On the other hand, the situation is subtle for Jinter- 
While Jinter is estimated to be 40 meV (> 0, i.e., AF) 
for structure A, Jinter ~ meV (paramagnetic) for B. 
This implies that F and AF interactions are both small 
in magnitude, so that the value of Jinter should depend 
on details of the crystal structure. 




FIG. 3. (a) The optimized atomic configuration of the 
crystallized polyaminotriazole doped with CIO4. The green, 
white, blue, red, pink, orange balls represent the position of 
the C, H, N, O, CI and F atoms, respectively, (b) The band 
structure for the ferromagnetic solution. The black (green) 
lines represent the bands for the majority (minority) spin, 
(c) The wave function (sum of the squared absolute values of 
the top-four, majority-spin valence bands at T) for the fer- 
romagnetic polyaminotriazole crystal doped with CIO4. Blue 
(red) contours represent the amplitude around PAT (dopant). 

This is why we move on to first-principles calculation of 
the flat-band ferromagnetism for actual chemical dopants 
introcuded in PAT crystal. We should opt for anions 
having large electron affinities, because the bands of the 
dopants to which the charge transfer from PAT occurs 
will then lie well away from the bands around the Fermi 
energy (i.e., the flat bands). In the present study, we 
have studied halogens (CI, F) and halogen compounds 
(CIO4, HF2). To make the flat bands half- filled, we put 
four dopants in each unit-cell that contains four five- 
membered rings. As for the polymorph we focus here 
on structure A, and we have adopted a slightly larger 
unit-cell size (9. 52 Ax 7. 41 Ax 7.17A) to accommodate 
CIO4 and HF 2 . 

We first found that we can exclude the Cl-doped case 
since the paramagnetic solution has a total energy lower 
than that for the ferromagnetic one by ~ 270 meV. For 
CIO4, F and HF2, on the other hand, the ferromagnetic 



state becomes lower in energy than the paramagnetic 
state. In Fig. 3 (a), we show the optimized atomic con- 
figuration for the ferromagnetic solution in the CIO4- 
doped crystal. We can see that the insertion of the CIO4 
molecules significantly changes the relative position of 
the two PAT chains in the unit-cell from those in Fig.l, 
which suggests that the molecules may crucially affect the 
electronic properties and unfavors the intra-chain ferro- 
magnetism. 

So let us turn to the band structure in Fig. 3(b) to see 
whether the states around Ep are modified. If the elec- 
tron affinity of CIO4 were sufficiently larger than that of 
PAT, the C104-charactered bands would not hybridize 
with the states around Ep. However, the result shows 
that the hybridization occurs for the flat (valence-top) 
bands, and the Mielke-Tasaki's orbits are severely af- 
fected. 
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FIG. 4. A plot similar to Fig. 3 for the fluorine-doped case. 

For a more quantitative argument, we have looked at 
the charge distribution by adopting a method due to 
Aizawa and Tsuneyuki [15,16] for estimating the Mul- 
liken charge for plane- wave basis band calculations. We 
divide the whole system into subsystems with each sub- 
system containing one atom, and calculate the valence 
charge for each. Summing the valence charges for dopant 
(CIO4 here), we can estimate the charge transfer from 
PAT to the dopant (Q). For PAT-CIO4, we have found 
that Q for each dopant molecule (i.e., for each five- 
membered ring) is c± 0.88, i.e., the charge transfer is in- 
complete. This suggests that the bands (Mielke-Tasaki's) 
to which we wanted to hole-dope are hybridized with the 
dopant (CIO4) band. We can indeed confirm this from 
the plot of the (majority-spin) wave function as a sum 
(p hereafter) of the squared absolute values of flat bands 
at T in Fig. 3(c). We can see that the amplitude extends 
substantially to CIO4, especially on O sites. In fact, the 
ground state becomes antiferromagnetic even within the 
chain, with Jintra = 60 meV and Jinter = 20 meV. 



Let us now move on to the fluorine-doped case. We 
show the optimized atomic configuration and the band 
structure in Fig. 4. The band structure is similar to those 
plotted in Fig. 2 with the charge transfer estimated to be 
Q ~ 0.78. The wave function in Fig. 4(c) indicates that 
the amplitudes on C atoms becomes considerably small, 
which means that the connectivity condition, necessary 
for the flat-band ferromagnetism [7], is violated. The ex- 
change energies are estimated to be Jmtra = 40 meV and 
Winter = 20 meV, so that the ground state is antiferro- 
magnetic both within and across the chains. 

Let us finally discuss the result for HF2, which has 
the largest electron affinity among the dopants studied 
here. While the band structure and wave function in 
Fig. 5 suggest that Mielke- Tasaki' s states hybridize with 
anion states to some extent, the charge transfer is cal- 
culated to be Q = 1.1, so that this case is more promis- 
ing than that of CIO4 or F. While Jmtra — 10 meV and 
Jintra — 30 meV are still both antiferromagnetic, Jmtra 
becomes considerably smaller than those for CIO4 and 
F, reflecting the high electron affinity of HF2. From the 
chemical trend, we may expect that a favorable situation 
(similar to the uniform doping described in the first half 
of the paper) for the ferromagnetism is expected to be 
realized for BF4 or PF6, which are known to have higher 
electron affinity, although a first-principles calculation for 
BF4 or PF6 would be too demanding, since these have 
many valence electrons or require large cut-off energies 
in the plane-wave expansion. 
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FIG. 5. A plot similar to Fig. 3 for the HF2-doped case. 

To summarize, a spin density functional calculation 
shows that the three-dimensional crystal of chains of five- 
membered rings has a robust intra-chain ferromagnetism 



originally conceived for an isolated polymer. For the 
chemical doping BF4 or PFg with high electron affin- 
ity should be promising. The magnetism in the organic 
crystal may also be controlled by the pressure effect. 
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